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1. Introduction 



There is no need to underline the success of the renormalization group revisited by 
Wilson p] [2] which is nowadays seen as a fundamental symmetry in lattice statistical 
mechanics or field theory. It contributed to promote 2-d conformal field theories 
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and/or scaling limits of second order phase transition in lattice statistical mechanicf0. 
If one does not take into account most of the subtleties of the renormalization 
group, the simplest sketch of the renormalization group corresponds to Migdal- 
Kadanoff decimation calculations where the new couplingconstants created at each 
step of the (real-space) decimation calculations are forcedj to stay in some (slightly 
arbitrary) finite dimensional parameter space. This drastic projection may be justified 
by the hope that the basin of attraction of the fixed points of the corresponding 
(renormalization) transformation in the parameter space is "large enough" . 

One heuristic example is always given because it is one of the very few examples 
of exact renormalization, the renormalization of the one-dimensional Ising model 
without a magnetic field. It is a straightforward undergraduate exercise to show 
that performing various decimations summing over every two, or three or ... N spins, 
one gets exact generators of the renormalization group reading Jjv : t — > t N where 
t is (with standard notations) the high temperature variable t = tanh(if). It is 
easy to see that these transformations Tjv, depending on the integer TV, commute 
together. Such an exact symmetry is associated with a covariance of the partition 
function per site Z{t) = C(t) ■ Z(t 2 ). In this particular case one recovers the (very 
simple) expression of the partition function per site, 2 cosh(if ), as an infinite product 
of the action of (for instance) T2 on the cofactor C(t). In this very simple case, this 
corresponds to using the identity (valid for |x| < 1): 



n 

n=0 
+3 



1 +x r 



1 



For T3 : t —> t one must use the identity 

n (!+*»" +*»»■) = n(V 



x 3 " +1 



n=0 n=0 



1 



and for Tjy : t — > t a similar identity where the 3 in the exponents is changed 
into N. 

Another simple heuristic example is the one-dimensional Ising model with a 
magnetic field. Straightforward calculations enable to get an infinite number of 
exact generators of the corresponding renormalization group, represented as rational 
transformation^: 

Tjy : (x, z) — > T N (x, z) = (x N , z N ) (1) 

where the first two transformations T2 and T3 read in terms of the two (low- 
temperature well-suited and fugacity-like) variables x — e 4 and z = e 2H : 

(x + z) (1 + xz) (1 + xz) 

x 2 = 7^— v> > Z 2 - ' 



and: 



(1 + z) 2 ' * x + z ' 

(z 2 x + 2z + l) (z 2 + 2z + x) _ z 2 x+2z+l 

{z 2 x + z + xz + x) z z + 2z +x 

% The renormalization group approach of important problems like first order phase transitions, 
commensurate-incommensurate phase transitions, or off-critical problems being more problematic, 
f In contrast with functional renormalization group 3, 4] \E\. 

I One simply verifies that these transformations reduce to the previous Tjv : t -> t N in the z = 1 
limit (no magnetic field). 
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One simply verifies that these rational transformations of two (complex) variables 
commute. This can be checked by formal calculations for Tn and Tm for any ./V 
and M less than 30, and one can easily verify a fundamental property expected for 
renormalization group generators: 

Tn ■ Tm = T M ■ T N — T NM , (2) 

where the "dot" denotes the composition of two transformations. The infinite number 
of these rational transformations of two (complex) variables (JT|) are thus a rational 
representation of the positive integers together with their product. Such rational 
transformations can be studied "per se" as discrete dynamical systems, the iteration of 
any of these various exact generators corresponding to an orbit of the renormalization 
group. 

Of course these two examples of exact representation of the renormalization group 
are extremely degenerate since they correspond to one-dimensional models^. Migdal- 
Kadanoff decimation will quite systematically yield raizona0 transformations similar 
to ([1]) in two, or more, variable4|f Consequently, they are never (except "academical" 
self-similar models) exact representations of the renormalization group. The purpose 
of this paper is to provide simple (but non trivial) examples of exact renormalization 
transformations that are not degenerate like the previous transformations on one- 
dimensional modelf@. In several papers [6j [7] for Yang-Baxter integrable models with 
a canonical genus one parametrization [8j [9] (elliptic functions of modulus k), we 
underlined that the exact generators of the renormalization group must necessarily 
identify with the various isogenies which amounts to multiplying or dividing r, the 
ratio of the two periods of the elliptic curves, by an integer. The simplest example is 
the Landen transformation [7] which corresponds to multiplying (or dividing because 
of the modular group symmetry r O 1/t), the ratio of the two periods: 

k ~> kL = "TXT' T 2t - ( 3 ) 

1 + k 

The other transformation^ correspond to t -h> N ■ t, for various integers N. In the 
(transcendental) variable r, it is clear that they satisfy relations like @. However, in 
the natural variables of the model (e K , tanh(iT), k = sinh^if), not transcendental 
variables like r), these transformations are algebraic transformations corresponding in 
fact to the fundamental modular curves. For instance ([3]) corresponds to the genus 
zero fundamental modular curve 

f ■ f ~ (3 + f)-(f + 1487 ■jj'+j' 2 ) 

+ 3 • 15 3 • (16 j 2 - 4027 jj' + 16 j' 2 ) (4) 
- 12 • 30 6 ■ {j + j') + 8 ■ 30 9 = 0, 

% For instance the fixed points of |[TJ are not isolated fixed points but lie on (an infinite number) of 
genus zero curves ... 

ff In well-suited Boltzmann weight variables like x and z in Q, and not in (bad) variables like K, 
the coupling constants or the temperature. 

|| Such representations of the renormalization group are not exact representations (the exact 
transformation acts in an infinite number of parameters) but some authors tried to define "improved" 
renormalization transformations imposing the compatibility (commutation) of the renormalization 
transformations with some known exact symmetries of the model (Kramers- Wannier duality, gauge 
symmetries ...). 

§ For which the partition function or other physical quantities are algebraic functions, 
f See for instance (2.18) in [lO]. 
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or: 

5 9 v 3 u 3 - 12-5 6 u 2 v 2 ■ (u + v) + 375 uv (16u 2 + 16v 2 - 4027vw) 
-64:(v + u) ■ (v 2 + US7vu +u 2 ) + 2 12 -3 3 -uu = 0, (5) 
which relates the two Hauptmoduls u = 12 3 /j(k), v — \2 3 / j(kif): 

(l-fc 2 + fc 4 ) 3 



j(k) = 256 



fc 4 • (1 - fc 2 ) 2 ' 

(l + 14fc 2 + fc 4 ) 3 



i(^)= is (1 _ fc2)4 . fc2 ■ 

One verifies easily that (j4]) is verified with j = j(k) and j' = j{ki,)- 

The selected values of fc, the modulus of elliptic functions, k = 0, 1 are actually 
/jxed points of the Landen transformations. The Kramers- Wannier duality k <-> 1/fc 
maps fc = onto fc = oo. For the Ising (resp. Baxter) model these selected values 
of fc correspond to the three selected subcases of the model (T — oo, T — and 
the critical temperature T = T c ), for which the elliptic parametrization of the model 
degenerates into a rational parametrization [7|. We have the same property for all 
the other algebraic modular curves corresponding to r ■<-> N ■ r. This is certainly 
the main property most physicists expect for an exact representation of a generator 
of the renormalization group, namely that it maps a generic point of the parameter 
space onto the critical manifold (fixed points). Modular transformations are, in fact, 
the only transformations to be compatible with all the other symmetries of the Ising 
(resp. Baxter) model like for instance, the gauge transformations, some extended 
s/(2) x sl(2) x sl(2) x sl(2) symmetry [TT], etc. It has also been underlined in [6, 7. that 
seeing ([3]) as a transformation on complex variables (instead of real variables) provides 
two other complex fixed points which actually correspond to complex multiplication 
for the elliptic curve, and are, actually, fundamental new singularities discovered on 
the x linear ODE [T31 [2J US] . In general, this underlines the deep relation between 
the renormalization group and the theory of elliptic curves in a deep sense, namely 
isogenics of elliptic curves, Hauptmodulqfl, modular curves and modular forms. 

Note that an algebraic transformation like ([3]) or ([5]) cannot be obtained from any 
local Migdal-Kadanoff transformation which naturally yields rational transformations: 
an exact renormalization group transformation like ([3]) can only be deduced from 
non-local decimations. The emergence of modular transformations as representations 
of exact generators of the renormalization group explains, in a quite subtle way, the 
difficult problem of how renormalization group transformations can be compatible with 
reversibility (iteration forward and backwards) . An algebraic modular transformation 
([S]) corresponds to r — > 2t and r — > r/2 in the same time, as a consequence of the 
modular group symmetry r <-> 1/r. 

| Suggesting an understanding 7 12 of the quite rich structure of infinite number of the singularities 
of the in the complex plane from a Hauptmodul approach [7] 1121 . Furthermore the notion of 

Heegner numbers is closely linked to the isogenies mentioned here [7] . An exact value of the j'-function 
j(r) corresponding one of the first Heegner number is, for instance, j(l + i) = 12 3 . 
f It should be recalled that the mirror symmetry found with Calabi-Yau manifolds 16 17 18 19 20 
can be seen as higher order generalizations of Hauptmoduls. We thus have already generalizations of 
this identification of the renormalization and modular structure when one is not restricted to elliptic 
curves anymore. 

X The fact that the renormalization group must be reversible has apparently been totally forgotten 
by most of the authors who just see a semi-group corresponding to forward iterations converging to 
the critical points (resp. manifolds). 
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A simple rational parametrizatior ljjl of the genus zero modular curve {S) reads: 

z z 2 /2 12 \ 

u = 1728 t ■ u = 1728 ~ = u( — . (6) 

(z + 16) 3 (z + 256) 3 V z J 

Note that the previously mentioned reversibility is also associated with the fact that 
the modular curve ([5]) is invariant by u <-> v, and, within the previous rational 
parametrization (|6]) , with the fact that permuting u and v correspond^ to the Atkin- 
Lehner involution z f-> 2 12 /z. 

For many Yang-Baxter integrable models of lattice statistical mechanics the 
physical quantities (partition function per site, correlation functions, ...) are solutions 
of selectecjjj] linear differential equations. For instance the partition function per site 
of the square (resp. triangular, etc.) Ising model is an integral of an elliptic integral 
of the third kind. It would be too complicated to show the precise covariance of 
these physical quantities with respect to (algebraic) modular transformations like ([5]). 
Instead, let us give, here, an illustration of the non-trivial action of the renormalization 
group on some elliptic function that actually occurs in the 2-D Ising model: a weight- 
one modular form. This modular form actually, and remarkably, emerged [22 in a 
second order linear differential operator factor denoted Zi occurring [13] for x (3) , and 
that the reader can think as a physical quantity solution of a particular linear ODE 
replacing the too complicated integral of an elliptic integral of the third kind. Let us 
consider the second order linear differential operator (D z denotes d/dz): 

, (z 2 + 56z + 1024) 240 
a = Di + -— -4z ■ D z 



(z + 16)(z + 64) z- (z + 16) 2 (z + 64)' 

which has the (modular form) solution: 

2 Fi([l/12,5/12],[l];1728 ^ + ' 16)3 ) (7) 

Do note that the two pull-backs in the arguments of the same hypergeometric function 
are actually related by the modular curve relation ([5]) (see ©). The covariance ([7]) 
is thus the very expression of a modular form property with respect to a modular 
transformation (t O 2t) corresponding to the modular transformation ([3]). 

The hypergeometric function at the rhs of ([7]) is solution of the second order linear 
differential operator 

. z 2 +416 z+ 16384 „ 60 
P = D 2 Z + - — • D z 



(z + 256) (z + 64) z z (z + 64) (z + 256) 2 ' 

which is the transformed of operator a by the Atkin-Lehner duality z «-» 2 12 /z, and, 
also, a conjugation of a: 

( z + 16 ^-1/4 ( z + 16 y/4 

ft Corresponding to Atkin-Lehnor polynomials and Weber's functions. 

§ Conversely, and more precisely, writing 1728 z 2 /(z + 256) 3 = 1728 z'/(z' + 16) 3 gives the Atkin- 
Lehner 21 involution z-z' = 2 12 , together with the quadratic relation z — z z' —48zz' — 4096 z' = 
0. 

1 They are not only Fuchsian, the corresponding linear differential operators are globally nilpotent 
or G-operators |22| . 
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Along this line we can also recall that the (modular form) functional: 

F(j)= r Vl2 • 2FJ [1/12, 5/12], [1];— V (9) 



verifies: 



z 

A relation like (J7J is a straight generalization of the covariance we had in the 
one-dimensional model Z(t) = C(t) ■ Z(t 2 ), which basically amounts to seeing 
the partition function per site as some "automorphic function" with respect to the 
renormalization group, the simple renormalization group transformation t —> t 2 being 
replaced by the algebraic modular transformation ([3]) corresponding to r O 2 r (that 
is the Landen transformation ©). 

We have here all the ingredients for seeing the identification of exact algebraic 
representations of the renormalization group with the modular curves structures we 
tried so many times to promote (preaching in the desert) in various papers [fU [?]■ 
However, even if there are no difficulties, just subtleties, these Ising-Baxter examples 
of exact algebraic representations of the renormalization group already require some 
serious knowledge of the modular curves, modular forms and Hauptmoduls in the 
theory of elliptic curves, mixed with the subtleties naturally associated with the various 
branches of such algebraic (multivalued) transformations. 

The purpose of this paper is to present another elliptic hypergeometric function 
and other much simpler (Gauss hypergeometric) second order linear differential 
operators covariant by infinite order rational transformations. 

The replacement of algebraic (modular) transformations by simple rational 
transformations will enable us to display a complete explicit description of an exact 
representation of the renormalization group that any graduate student can completely 
dominate. 

2. Infinite number of rational symmetries on a Gauss hypergeometric 
ODE 

Keeping in mind modular form expressions like (0, let us recall a particular Gauss 
hypergeometric function introduced by R. Vidunas in |25) : 

This hypergeometric function corresponds to the integral of a holomorphic form on a 
genus- one curve P(y, t) = 0: 

-, with: y 4 - £ 3 ■ (1 - £) 2 = 0. (12) 

y 

Note that the function 

Hz) = z lfi - 2 F 1 ([ii],[J];z), (13) 

| Where j is typically the j'-function 1101 123) . 
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which is exactly an integral of an algebraic function, has an extremely simple 
covariance property with respect to the infinite order rational transformation z 
-4z/(l-z) 2 : 

^TT%) = (-4) 1/4 -^)- (14) 



-(1 

The occurrence of this specific infinite order transformation is reminiscent of Kummer's 
quadratic relation 

2 F 1 ([a, b],[l + a-b}; z) = (15) 
= (l-z)- a - 2 Fi ([-,-I__6] > [l + o-6]; 



V L 2' 2 ~ J ' L ~ ' " {l-z)\ 

but it is crucial to note that, relation (fT4|) does not relate two different functions, but 
is an "automorphy" relation on the same function. 

It is clear from the previous paragraph that we want to see such functions as 
'ideal' examples of physical functions covariant by an exact (here, rational) generator 
of the renormalization group. The function (|13[) is actually solution of the second 
order linear differential operator: 

°= D l + 7 3 n 5 —\ ■ D * = with: ( 16 ) 

4 z ■ (1 — z) 

1 3-5z „ 1 dln(z 3 (1 - z) 2 ) 

«x = D„ + - — r = D z 1 



4 z ■ (1 - z) 4 dz 

From the previous expression of uj\ involving a log-derivative of a rational function 
it is obvious that this second order linear differential operator has two solutions, the 
constant function and an integral of an algebraic function. Since these two solutions 
behave very simply under the infinite order rational transformation z — > — 4z/(l— z) 2 , 
it is totally and utterly natural to see how the linear differential operator f2 transforms 
under the rational change of variable z — > R{z) = — 4z/(l — z) 2 (which amounts to 
seeing how the two order-one operators oj\ and D z transform) . It is a straightforward 
calculation to see that introducing the cofactor C(z) which is the inverse of the 
derivative of R{z) 

C[Z) - 4 l + z ' C{z) ~ dz ' (17) 

D z and w\ respectively transform under the rational change of variable z — > R{z) = 
-4z/(l - z) 2 as: 

D z — > ^(z).^, wi — > = C{zf ■u Jl --^—, (18) 

o(_zj 

yielding: fi — ► C*(z) 2 • fi. (19) 

Since z — > — 4z/(l — z) 2 is of infinite order, the second order linear differential 
operator (|16[) has an infinite number of rational symmetries (isogenics) : 

- Az i 16. (1 " Z)2 ' Z 



(1-z) 2 (1 + ^) 4 

( l-z) 2 (l + z) 4 
(1 -6z +z 2 ) 4 

Once we have found a second order linear differential operator (written in a 
unitary or monic form) f2, covariant by the infinite order rational transformation 
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z — > — 4z/(l — z) 2 , it is natural to seek for higher order linear differential operators 
also covariant by z —> — iz/(l — z) 2 . One easily verifies that the successive symmetric 
powers of f2 are (of course ...) also covariant. The symmetric square of fi, 

3 3 — 5 z 2 3 1 — 5 z 

4 (1 - z) z z + 8 (l-z)z 2 
factorizes in simple order-one operators: 

2 3-5z \ / 1 3-5z 



D l +l^—^- D l +Z n W -D„ (21) 



4 (1 - z)z/ V 4 (1 - z) z 

and, more generally, the symmetric iV-th powejj] of Q reads 

N 3-5z \ / JV-1 3-5z \ 

4 z(l-z)/ \ 4 z(l - z)J 

1 3-5z\ n 
4 z(l — z)J 

The covariance of such expressions is the straight consequence of the fact that the 
order-one factors 

= D z + - 3 ; 5Z fc = 0, 1, N, (24) 

transform very simply under z — > —Az/(1 — z) 2 : 

0Jk -> K) W = (C(z)) fc+1 - Wfe .(c(z))^. (25) 

More generally, let us consider a rational transformation z — > R(z), the corresponding 
cofactor C(z) — 1/R'(z), and the order-one operator ui = D 2 + We have 

the identity: 

' ^ ' %) J * dT - * (26) 

The change of variable z — > R(z) on oj-y reads: 

D z + A(z) — ► C(«).£>, + A(i?(z)) - C{z) ■ (D, + B(z)). 
We want to impose that this rhs expression can be written (see l[TS)l ) as: 

C(z) 2 • (p. + A(z)) ■ 

which, because of ([26)) . occurs if 

B(z)= A(z)- dl ^\ 
yielding a "Rota-Baxter-like" [2l>ll2"7] functional equation on A(z) and R(z): 

(^Y • A{R{z)) = d -^l-A(z) +^#. (27) 
V dz ) dz dz 

Remark: Coming back to the initial Gauss hypergeometric differential operator 
the covariance of Q becomes a conjugation. Let us start with the Gauss 
hypergeometric differential operator for ([TTj) : 

H = 8z- (1-z) -D 2 Z + 2- (5 -7z) ■ D z - 1. (28) 

f Such formula is actually valid for Qa = (Dz ' D z for any A(z). Denoting 5jv symmetric 

JV-th power of Qa one nas <Sjv = (-Dz + A(z)) ■ cSjv-i- 
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It is transformed by z —> R(z) — — 4z/(l — z) 2 into: 

Hi = 8z ■ (1 - z) ■ D 2 Z - 2 (3 z - 5) • D z + -— — 

= (l-^.tfo-CL-*)- 1 ^, 
then by z -> R{R{z)) = R 2 {z) = 16 z (1 - z) 2 /(l + z) A into: 

, f3z-l)(z + 5) n , 1C z-1 



iJ 2 = 8 z ■ (1 - z) • Di - 2 - '-- -D z +16 



z + 1 

z + 1 \ TT / z + 1 \-l 

#0 



(^ + 1) 2 



and more generally for z — > ii/v = R(R{R ■ ■ ■ {R(z) ■ ■ ■): 

H N = Cn-Hq-C^ 1 , where: C N = z 1/4 • i?^ 74 . 

1. A few remarks on the "Rota-Baxter-like" functional equation. 

The functional equation^] (f2"T)l is the (necessary and sufficient) condition for fl 
(D z + A(z)) ■ D z to be covariant by z — > R(z). 

Using the chain rule formula of derivatives of composed functions: 

dR(R(z)) dR(z) .dR{z) 



dz dz dz 

d 2 R(R(z)) 



-(R(z))} 



dz 2 

d 2 R(z) ,dR{z) (dR{z)\i d 2 R(z 



(*(*))]+ (-^) 



dz 2 dz 

one can show that, for A(z) fixed, the "Rota-Baxter-like" functional equation (j2"T)l is 
invariant by the composition of R(z) by itself -R(z) — > R(R(z)), R(R(R(z))), ■ ■ ■ 
This result can be generalized to any composition of various R(z)'s satisfying (|27|) . 
This is in agreement with the fact that (|27|) is the condition for $1 = [D z + A(z)) ■ D z 
to be covariant by z — > R(z): it must be invariant by composition of i?(z)'s (for A(z) 
fixed). 

Note that we have not used here the fact that for globally nilpotent 22 operators, 
A(z) and B(z) are necessarily log-derivatives of N-th. roots of rational functions. For 
R(z) = -4z/(l-z) 2 : 

= 1 dHa(z))^ = I dHHz)) 

A dz A dz 

a(z) = (l-z) 2 -z 3 , b(z) = z 3 {1 + Z)4 



(1-z) 10 ' 

The existence of the underlying a(z) in (|2U)) consequence of a global nilpotence of 
the order-one differential operator, can however be seen in the following remark on 
the zeros of the lhs and rhs terms in the functional equation (|27|) . When R{z) is a 
rational function (for instance — 4z/(l — z) 2 or any of its iterates R^(z)), the lhs 
and rhs of (f27|) are rational expressions. The zeros are roots of the numerators of these 



f The Rota-Baxter relation of weight © reads: R(x) R(y) + 0R(x y) = R(R(x) y + x R(y)). 
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rational expressions. Because of (|2T))) the functional equation (|2"T)) can be rewritten 
(after dividing by R'(z)) as: 

. + ,30, 

One easily verifies, in our example, that the zeros of the rhs of (|30|) come from the 
zeros of A{R(z)) (and not from the zeros of R'(z) in the lhs of (|30p). The zeros of 
the log-derivative rhs of (|30|) correspond to a(z) • R'(z) 4 — p, where p is a constant 
to be found. Let us consider for R(z) the n-th iterates of — 4z/(l — z) 2 that we 
denote R^(z). A straightforward calculation shows that the zeros of A(R^(z)) or 
a'(R( n )(z)) (where a'(z) denotes the derivative of a(z) namely (z — 1) (5 z — 3) • z 2 ) 
actually correspond to the general closed formula: 

' 5 ^ z) ' ( dRi 2 {z) y ~ 4 - 3 M- 4 )" = °- ( 3i ) 



5 ;i 



More precisely the zeros of 5-R^(z)-3 verify (|3Tj), or, in other words, the numerator 
of 5 R^ n \z) — 3 divides the numerator of the lhs of (|3"Tj) . 

In another case for T(z) given by (|43| . which also verifies ([27]) (see below), the 
relation ([31]) is replaced by: 

55 - Q ( z )'( dz ) " 4-3 3 -(-7 - 24i) n = 0. (32) 

More generally for a rational function /o(x), obtained by an arbitrary composition of 
— 4z/(l — z) 2 and T(z), we would have: 

5 5 . ffl (z)- (^) 4 ~ 4-3 3 -A" = 0. (33) 
where A corresponds to: 

p(x) = A • z + ■ ■ ■ , A = . (34) 

L dz J 2=0 

Symmetries of 57, solutions the "Rota-Baxter-like" functional equation. 

Let us now analyse all the symmetries of the linear differential operator 51 = 
(D z + A(z)) ■ D z by analyzing all the solutions of (|27|) for a given A(z). For 
simplicity we will restrict to A(z) = (3 — 5z)/z/(l — z)/A which corresponds to 
R(z) = —4z/(z — l) 2 and all its iterates (|2"0"1) . Let us first seek for other (more 
general) solutions that are analytic at z = 0: 

R{z) = ai ■ z + a2 ■ z 2 + <23 • z 3 + • • • (35) 

It is a straightforward calculation to get, order by order from (|2"7) . the successive 
coefficients a n in (|35|) as polynomial expressions (with rational coefficients) of the 
first coefficient a\ with 

2 1 
a2 = -r • oi • (ai - 1), a 3 = — • ai • (oi - 1) • (7ai - 17), 

2 

"4 = -7^7 ' a i ' Oi - !) ' ( 41 a i - 232 fl i + 366 )> ' ' ' 
4875 

a n = - - ■ oi ■ (ai - 1) • 5 ; tt , (36) 
5 -P«(-4) 
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R 



where P n (ai) is a polynomial with integer coefficients of degree n — 2. Since we have 
here a series depending on one parameter a\ we will denote it R ai (z). This is a 
quite remarkable series depending on one parameteiQ- One can easily verify that this 
series actually reduces (as it should!) to the successive iterates (|2"U1) of— 4z/(l — z) 2 
for a\ — (—4)™. In other words this one-parameter family of "functions" actually 
reduces to rational functions for an infinite number of integer values a\ — (—4)™. 

Furthermore, one can also verify a quite essential property we expect for a 
representation of the renormalization group namely that two R ai (z) for different 
values of a\ commute, the result corresponding to the product of these two a\\ 

Ql (i? hl (z)) = R bl (R ai {z)) = R ai - bl {z). (37) 

The neutral element must necessarily correspond to a\ = 1 which is actually the 
identity transformation R\(z) = z. We have an "absorbing" element corresponding 
to a\ = 0, namely Ro{z) = 0. Performing the inverse of R ai (z) (with respect to 
the composition of functions) amounts to changing a\ into its inverse l/a±. Let us 
explore some "reversibility" property of our exact representation of a renormalization 
group with the inverse of the rational transformations (f2"U| . The inverse of R-4 (z) = 
— 4z/(l — z) 2 must correspond to a\ = —1/4: 

*-•/.<»>-- l-t-r 2 -^-^-^ + - < 38 > 

However, a straight calculation of the inverse of R- 4 (z) = — 4z/(l — z) 2 gives a 
multivalued function, or if one prefers, two functions: 

*- 2+ l^ = -i.. (39) 



^(21 z-2 -2y/l- Z 4 „ 1 1 o 

S% 4 (z) - ^ = _ +2 +-z +-z 2 + • • • , 

1/4 z z 4 8 

which are the two roots of the simple quadratic relation (i?_4(z') = z): 

z 12 - 2-(l -~W +1 = 0, (40) 

z 

where it is clear that the product of these two functions is equal to +1. The radius 
of convergence of S_y 4 (z) is R = 1. 

Because of our choice to seek for functions analytical at z = our renormalization 
group representation "chooses" the unique root that is analytical at z = 0, namely 
S^l, 4 (z). For the next iterate of R-i(z) = -4z/(l — z) 2 in (|20| the inverse 
transformation corresponds to the roots of the polynomial equation of degree four 
(Rie(z') = z): 

z * +^--)-z' 3 +(6 + -)-z' 2 +(4- -)-z' +1 = 0, (41) 

z z z 

which yields four roots, one of which is analytical at z = and corresponds to 
ai = l/(— 4) 2 in our one-parameter family of (renormalization) transformations: 

(1) _L_ J{_ 2 ^3_ 3 277 4 3181 5 

■^i/iei 2 ) - i 6 z + 128 Z + 4096 2 + 32768 Z + 524288 Z + ' ' ' 

% For A(z) given we get a one-parameter family of R(z) solution of I I27II , Conversely, for R(z) given 
one can ask if there are several A(z) such that H27t is verified. This is sketched in |Appcndix A| 
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its (multiplicative) inverse Sy 16 (z) — 1/ S^ le (z): 

q( 2 ) r ^ - 16 _ x 17 7445 4 

^VieW - z b i 6 z i28 Z 4096 Z 32768 Z + "' 

and two (formal) Puiseux series (u — ±y/z): 

S\ 3 Uz) = 1 +u +-u 2 +-u 3 +-u i + —u 5 + —u e + •■• 
1 /!6 V ; 2 8 4 128 32 

Many of these results are better understood when one keeps in mind that there is a 
special transformation J : z 1/z which is also a R-solution of (|27p and verifies 
many compatibility relations with these transformations (Id denotes the identity 
transformation Rq(z)): 



R—4 • J — -/?_4, ^ 1/4 ' -^—4 — '^i 



c(l) - <?W c(2) / >, _. c(l) o(2) 

''l/ieW — D -l/4 ' J -l/4' D l/ieW — D -l/4 ' D -l/4> 

T o(2) t c(2) _ o(l) 

J ' J -l/4 — °-l/4' J ' J -l/4 — °-l/4> 

where the dot corresponds, here, to the composition of functions. These symmetries 
of the linear differential operator f2 correspond to isogenies of the elliptic curve (fT2")l . 

It is clear that we have another one-parameter family corresponding to J ■ R ai 
with an expansion of the form: 

2 (h -l)(4 6i +l)(4 6i +3) 2 



975 6 2 



2 



1 (6i - 1) (4 h + 1) (1268 b\ + 951 &i + 91 



z 



248625 fef 

2 (&! - 1) (4 h + 1) (3688 6? + 2766 b\ + 404 b x + 17) 
~~ 2071875 b\ 
H 

For b a = -1/4, bi = (-1/4) 2 , b x = (-1/4) 3 , this family reduces to the 
(multiplicative) inverse of the successive rational functions displayed in (|2"0|) 

1 (1-z) 2 _t_ (1 + z) 4 

' 4 z ^ 16 ' (1 - z) 2 ■ z 

_L_ (1 -6z + z 2 ) 4 

~64 ' (l-z) 2 -(l + z) 4 -z ~ ^ '"' 

which can also be written as: 

1 . 1, 1 1 . 1. 3 z 

(* + -) + ^ T7'( z + -)+o + 



4 v z ; 2' 16 v z y 8 (1-z) 2 ' 

1 . 1, 13 z 17-60z+102z 2 -60z 3 + 17z 4 



64 v z y 32 4 (1-z) 2 (1 + z) 4 

1 I 51 ; z 17-60z+102z 2 -60z 3 + 17z 4 

256' (Z+ z' 1 + 128 + 16 (1-z) 2 (1 + z) 4 
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z- (1 -z) 2 (l + z) 4 



16 



(z 2 - 6z + l) 4 



1 , l s 205 z 17-60z + 102z 2 -60z 3 + 17z 4 
(z + -) + 



1024 v z' 512 164 (1 - z) 2 (1 + z) 4 

1 z-(l-z) 2 (l + Z ) 4 _ g4 z-(l-z) 2 (l + z) 4 (z 2 -6z+l)' 



(z 2 -6z+l) 4 (1 + 20 z - 26 z 2 + 20 z 3 + z 4 



' (^ + -) + ( 4 " 



(-4)™ v z ; 54™ 

z 17-60z+102z 2 -60z 3 + 17z 4 



(-4)™- 2 (l-z) 2 (l + z) 4 

z (l-z) 2 (l +z) 4 

+ (-4)"- 6 ' (z 2 -6z + l) 4 

z (l-z) 2 (l + z) 4 (z 2 -6z + l) 4 

+ (_4)«-8 ' (i + 20 z - 26 z 2 + 20 z 3 + z 4 ) 4 + '"' 

where we discover some "additive structure" of these successive rational functions. 

In fact, due to the specificity of this elliptic curve (occurrence of complex 
multiplication), we have another remarkable rational transformation solution of (|27[) . 
preserving covariantly Q. Let us introduce the rational transformation ( i denotes 

we also have the remarkable covariance [25] : 

/ r l 1, r 5, \ 1 - z/(l + 2i) (A l n .5, m/ . \ 

which can be rewritten in a simpler way on (|13l) (see (|14p ). 

It is a straightforward matter to see that T{z) actually belongs to the R ai (z) 
one-parameter family: 

T(z) = R ai (z) = -(7 + 24i)-z + oi = -25 -p, 

p = (7 + 24z)/25, \p\ = 1. 

As far as the reduction of (|36[) to a rational function is concerned, it is 
straightforward to see that: 

(l-z) 2 -(l + z) 4 -i? Ql (z) = 0LX-Z+--- (44) 
2 ■ oi • (oi - 1) • (oi + 4) • (oi - 16) • P 8 (ai) • z 8 



(42) 



175746796875 

+ ... 

- 777^ • oi • (oi - 1) • (oi + 4) • (ai - 16) • P n (ai) ■ z n 
iV (n) 

H 

where iV(n) is a large integer growing with n, and P„ is a polynomial with integer 
coefficients of degree n — 4, or 

(1 - (1 +2i)-z) 4 -i? Ql (z) = ai •*+••■ 
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- 124 3i25 ' ai ' (ai _1 )'( ai + 7 +24*)-(- p e(oi) -MQ 6 (ai))-^ 6 
+ • • ■ (45) 

+ TFTT ' a i ' (ai - 1) • («i +7 + 24i) • (P„(oi) + iQ n {a 1 ))-z n 
N(n) 

H 

where P n and <5 n are two polynomials with integer coefficients of degree respectively 
n — 3 and n — 4. 

Similar calculations can be performed for T*(z) defined by 

for which we also have the covariance: 

It is a simple calculation to check that any iterate of T(z), (resp. T*(z)) is 
actually a solution of (127[) and corresponds to P ai (z) for the infinite number of values 
a\ = (—7 — 24 i) N (resp. (—7 + 24 i) ). Furthermore, one verifies, as it should 
(see (|37)l). that the three rational functions i?-4(z), T(z), and T*(z) commute. It 
is also a straightforward calculation to see that the rational function built from any 
composition of R-^{z), T(z) and T*(z) is actually a solution of (|27[). We thus have 
a iripZe infinity of values of ai, namely ai = (— 4) M • (—7 — 24i) w • (—7 + 24 i) p for 
any integer M , N and P, for which R ai (z) reduces to rational functions. We are in 
fact describing (some subset of) the isogenies of the elliptic curve (fT2"|) . and identifying 
these isogenies with a discrete subset of the renormalization group. Conversely, a 
functional equation like (|27[) can be seen as a way to extend the n-fold composition 
of a rational function R{z) (namely R(R(- ■ ■ R(z) ■ ■ •))) to n any complex number. 

2.3. Revisiting the one-parameter family of solutions of the "Rota- Baxter-like" 
functional equation. 

This extension can be revisited as follows. Keeping in mind the well-known example of 
the parametrization of the standard map z — > 4z- (1 — z) with z = sin 2 (#), yielding 
9^26, let us seek for a (transcendental) parametrization z — P(u) such that 

R-Jp(u)\ = P(-iu) or: i?_ 4 = P • ■ P~\ (47) 

where H ai denotes the scaling transformation z — > a% ■ z (here 7J_4 : z — > —4 • z) 
and P _1 denotes the inverse transformation of P (for the composition). One can 
easily find such a (transcendental) parametrization order by order: 

2 2 7 3 82 4 1078 5 452 6 

(Z) ~ Z ~5 Z + 75 Z ~ 4875 Z + 414375 Z ~ 1243125 Z 

57311 7 1023946 H 

H z 7 z 8 + ■■■ 48) 

1212046875 175746796875 1 ; 

and similarly for its inverse (for the composition) transformation: 

n-i/ n 2 2 17 3 244 4 45043 5 

Q(*) + ~5 Z + 1625 Z + 414375 Z 

2302 6 128941 7 15365176 8 



Z + mono7t; Z + ooiin«7i; Z + (^) 



27625 1939275 281194875 
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This approach is reminiscent of the conjugation introduced in Siegel's theorem [281 
[221 [3D]- It is a straightforward matter to see (order by order) that one actually has 

R ai (P{u)) = P{ax-u) or: R ai = P ■ H ai ■ P' 1 . (50) 

The structure of the (one-parameter) renormalization group and the extension of the 
composition of n times a rational function R{z) (namely R(R(- ■ ■ R(z) ■ ■ ■))) to n any 
complex number, becomes a straight consequence of this relation. Along this line one 
can define some "infinitesimal composition" (e ~ 0): 

R 1 + £ (z) = p.H 1 + e -P-\z) = z + e-F(z) + (51) 

where one can find, order by order, the "infinitesimal composition" function F(z): 

r.r S 2 2 2 3 14 4 154 5 22 6 

F(Z)= -15* 3 "195* --mT - 663* 

418 - 7 -.!S^ 8 -4^ 9 +-" (52) 



16575 480675 160225 
It is straightforward to see, from (|3"7)l . that the function F(z) satisfies the following 
functional equations involving a rational function R(z) (in the one-parameter family 

RaM)- 

^■Fiz) = F(R(z)), d -^M.F(z) = F(R^z)), 
az az 

where: R {n \z) = R(R(- ■ ■ R(z)) ■ ■ •)• (53) 

F(z) cannot be a rational or algebraic function. Let us consider the fixed points of 

R( n \z). Generically is not equal to or oo at any of these fixed points. 

Therefore one must have F(z) — or F(z) = oo for the infinite set of these fixed 
points: F(z) cannot be a rational or algebraic function, it is a transcendental function, 
and similarly for the parametrization function P{z). In fact, let us introduce the 
function 

G(z) = (l-z)-F(z), (54) 

^, s 7 2 4 3 4 4 28 5 44 6 44 7 

G(z) = z z 2 H z H z 4 H z 5 H z° -\ z 7 

y ' 5 15 65 1105 3315 5525 

836 8 1748 n 

■ z + z 9 + ■■■ + g n ■ z n + • • ■ 



160225 480675 

One actually finds that the successive g n satisfies the very simple (hypergeometric 
function) relation: 

£n+l = 4n_^ 
g n An + 1 

The function G(z) is actually the hypergeometric function solution of the 
homogeneous operator 



Dl + ~ • D z 



1 13z - 3 „ 36z^-3z + l 



4 z-(l-z) 4 (l-z) 2 -z 2 ' 

or of the inhomogeneous ODE: 

4z-(l-z)-^ + (9z -3) -G(z) -z- (1-zf = 0. 
az 

One deduces the expression of F(z) as a hypergeometric function 

F(z)= z-(l- Z )^. 2 ^([i,i],[J;,) = ^k = i- (56) 
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Finally we get the linear differential operator annihilating F(z) 



tt F = £>?+-• — - D, 



1 5z-3 13-6z + 5 2 2 



4 z(l-z) 4 (l-z) 2 z 2 

= D.-fD.-l.^"), (57) 



4 z-(l-z) 

which is, in fact, nothing but fi* the adjoint of linear differential operator f2 (see 
(|16[) ). One easily checktlg that the second order differential equation Hf (y(z)) = 
transforms under the change of variable z — > — 4z/(l — z) 2 into the second order 
differential equation f2^(t/(z)) = with — C(z) 2 ■ ojp where the unitary 

(monic) operator uj f '' is the conjugate of dp- 

n 2 1 ll^ 2 + 30z + 3 

Up - Dz ~ r z-(i-z){i+z)' Dz (58) 

1 3 + 12z + 50z 2 + 12z 3 +3z 4 

+ 4 z 2 ■ (1 - z) 2 (1 + z) 2 

1 \ / 1 3-5z \ 

D Z .[D Z -- ■ ———).C(z) 



0. 



C(z)/ z V 2 4 z-(l-z), 

= (^))-^-^)= 

with C(z) = 1/R'(z) and the 'dot' denotes the composition of operators. Actually, 
the factors in the adjoint f2* transform under the change of variable z — > — 4 2/(1— z) 2 
as follow^: 

A* — > <?(*)• £> z , wf = — > = wf • 

ft* — > n^ fl) = C(z) • Q* • C(z) (59) 

which is precisely the transformation we need to match with (|53l) and see the ODE 
fl*(F(z)) = compatible with the change of variable z — > — 4z/(l — z) 2 : 

n* (F(z)) = — > (C{z) ■ n* ■ C{z)) (F(R(z))) (60) 

= (C(z) ■ ft* ■ C(zj)(B!(z) ■ F(z)) = C(z) ■ ft*(F(z)) 

This is, in fact, a quite general result that will be seen to be valid in a more general 
(higher genus) framework (see (I130|) . (I132|) below ). 

Not surprisingly one can deduce from (|37p and the previous results, in particular 
(|5"6"j) . the following results for R ai (z): 

-4.^| 01 = _ 4 = F(i?(z)), 

r) R 

(-4) n -^| ai = ( -4r = F(i2<">(*)), 

where R{z) = -4z/(l-z) 2 and R^{z) denotes R(R(- ■ ■ R(R(z))j). Of course we 
have similar relation for T(z), —4 being replaced by —7 — 24 i. Therefore the partial 

§ Using the command "dchange" with PDEtools in Maple. 

f Note that the result for lj* is nothing but transformation l|25ll on u)^ for k = — 1. Also note 
that the two transformations, performing the change of variable z — > — 4z/(l — z) 2 and taking the 
adjoint, do note commute: (ui^ y ^ ^ (iuj\y- R > Y ■ 
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derivative dR ai /dai that can be expressed in terms of hypergeometric functions for 
for a double infinity of values of ax, namely a\ = (— 4) M x (—7 — 24)^. 

One can, of course, check, order by order, that (|53p is actually verified for any 
function in the one-parameter family R ai {z): 

^^l.F(z) = F(R ai (z)). (61) 

which corresponds to an infinitesimal version of (|37[) . 

From (|5ip one simply deduces 

z . d l^± = F{ p {z% ( 62 ) 
dz 

that we can check order by order from (|48l) . the series expansion of P(z), and from 
(|52p the series expansion of F(z), but also 

^M-F(z) = Q(z) (63) 

that we can, check order by order, from (|49| . the series expansion of Q(z) = P^ 1 (z) 
and from (|52|). We now deduce that the log-derivative of the "well-suited change of 
variable" Q(z) is nothing but the (multiplicative) inverse of a hypergeometric function 

(fife) 

The function Q(z) is solution of the non-linear differential equation 

- 4 z 2 ■ (1 - z) 2 • (Q • Q (1) • Q< 3 > + (Q«) 2 • Q< 2 ) -2Q- (Q< 2 >) 2 ) 

+ z-(3-5z)(l-z)-Q«- (Q-Q^ - (Q (1) ) 2 ) 

+ (5z 2 -6z +3) Q- (Q (1) ) 2 = 0, (65) 

where the Q^'s denote the n-th derivative of Q(z). At first sight Q(z) would be a 
non-holonomic function, however, remarkably, it is a holonomic function solution of 
an order-five operator which factorizes as follows: 

^ 3-5z \ /■ 3 3-5z \ ,„„. 

(^ + (i-^)-(^ + 4'(T^)- (66) 

/ 2 3-5z \ /_ 1 3-5z \ _ 



(l-z)-z/ V 4 (l-z)-z 
yielding the exact expression of Q(z) in terms of hypergeometric functions: 

Q(z) = z- [F x {[\, I], [|; z)) 4 = (67) 

1 — z \2 44 4 1 — z / 
that is the fourth power of (TTBl . the differential operator ([66]) being the symmetric 
fourth power of f2. From (TIB")) we immediately get the covariance of Q(z): 
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and, more generally, Q^i? ai ^ = a± ■ Q(z). Since Q(z) and F(z) are expressed in 

term of the same hypergeometric function, the relation (|63[) must be an identity on 
that hypergeometric function. This is actually the case. This hypergeometric function 
verifies the ingomogeneous equation: 

4z dn^z} _(i_ z )-i/2 = o, (69) 

dz 

where: U(z) = 2 F 1 ([i i], [|; 

Recalling Q(P(z)) — z, one has the following functional relation on P(z) 

P(z)- 2Fi([\, i],[J];F(z)) 4 = z. (70) 

Noting that Q(z 4 ) 1 / 4 = F(z A ) (see (USD) can be expressed in term of an 
incomplete elliptic integral of the first kind of argument y— 1 

f J 1 ' Jl z 4 ) = EllipticF{z, V-i), (71) 



z ■ 2 Fi(\ ■ ■ 

V4 2 J L 4 J 

one can find that (170)1 rewrites on -P(z) as 

EllipticF{P{z) 1 / A , V^l) = z 1/A , (72) 
from which we deduce that the function P{z) is nothing but a Jacobi elliptic functio^ 

P{z) = ( S n(zy\ y^I)) 4 . (73) 

In |Appcndix B| we display a set of "Painleve-like" ODEiJj] verified by P(z). From the 

simple non-linear ODE on the Jacobi elliptic sinus, namely S"' -I- 2 • S* 3 = 0, and 

the exact expression of P{z) in term of Jacobi elliptic sinus, one can deduce other 

non-linear ODEs verified by the non-holonomic function P(z) (P^ = dP(z)/dz, 
P {2) = d 2 P ( z y dz 2y. 

^•(PW) 2 - (l-P)-P 3 / 2 = 0, (74) 

l2) 3 (PW) 2 3 i p3/2 

pi 2 ) . 2 £_ -l _ . l _ . = n (75) 

4 P + 4 z 2 2 3/ 2 1 j 

Singularities of the Jacobi elliptic function P(z). 

Most of the results of this section, and to some extent, of the next one, are straight 
consequences of the exact closed expression of P(z) in term of an elliptic function. 
Following the pedagogical approach of this paper we will rather follow a heuristic 
approach not taking into account the exact result (|73| , to display simple methods and 
ideas that can be used beyond exact results on a specific example. 

From a diff-Pade analyzis of the series expansion of P(z), we got the sixty 
(closest to z = 0) singularities. In particular we got that P(z) has a radius of 
convergence R ~ 11.81704500807 •• • corresponding to the following (closest to 
z = 0) singularity z = z s of P{z): 

z s = -11.817045008077115768316337283432582087420697 ••• 



= (-4)- ^([i |u|];i) 4 = —■ f £&- m 



16 r(3/4) 5 

% Denoted JacobiSN in maple: P(z) = (J acobiS N '{z 1 ' ' 4 , . 

f As a (non-holonomic) elliptic function P(z) provides elementary examples 1241 of non-linear ODEs 
with the Painleve property (like the Wcierstrass-P function). 
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This singularity corresponds to a pole of order four: P(z) ~ (z — z s )~ 4 . The function 
P(z) has many other singularities: 

3 4 -z s , (161 ± 240i) • z s , (-7 ± 24 i) • z s , (-119 ± 120 i) ■ z a , • 
5 4 -z s , (41 ± 84(H) • 2 a , (-527± 336i)-z s , (-1519 ± 720 i) • z s , 
7 4 -z s , (1241 ± 2520 i) ■ z s , (-567 ± 1944 i) ■ z s , (-3479 ± 1320 i) ■ z s 

In fact, introducing x and y the real and imaginary part of these singularities in 
z s units, one finds out that they correspond to the double infinity of points 

x = (m\ — 2 mi m 2 — m 2 ) • (to? + 2 mi TO2 — m 2 ), 

y = 4 • mi W2 • (m 2 — mi) ■ (m 2 + mi), (77) 

where mi and m 2 are two integers, and they all lie on the intersection of an infinite 
number of genus zero curves indexed by the fourth power of an integer M = m 4 
(m = mi or m = m 2 ): 

2 12 M 4 ~2 u -x-M 3 — 2 7 • (17 y 2 + lAx 2 ) ■ M 2 

- 2 5 -x- (8x 2 +7y 2 ) ■ M + y 4 = 0. (78) 

The parametrization (|77|) describes not only the poles of P(z) when mi + m 2 is odd, 
but also the zeros of P(z) when mi + m 2 is even. This (infinite) proliferation of 
singularities confirms the non-holonomic character of P(z). 

These results are simply inherited from ff73|) . The zeros and poles of the elliptic 
sinus sn(z, i) correspond to two lattice of periods. Denoting /Ci and /C 2 the two 
periods of the elliptic curve, the location of the poles and zeros reads respectively: 

P nu n 2 = 2n 1 -K 1 + (2n 2 + 1) ■ JC 2 , (79) 

Z ni ,n 2 = 2m ■ fCi + 2n 2 ■ /C 2 , (80) 

Kl = ■ ' = (1 - V=I) • *i (81) 

making crystal clear the fact that we have complex multiplication for this elliptic 
curve. The formula (|77p just amount to saying that the poles and zeros of sn(z 1 / 4 , i) 
are located at P 4 i n2 and Z 4 i n2 : 

P 4 

ni,n 2 



7 4 — —El 

n 1 ,n 2 ~ | 



j • ((2 m + 2 n 2 + 1) + i ■ (2 n 2 + 1)) 

((2m + 2n 2 + 1) + i ■ 2n 2 y. (82) 



The correspondence with (|77[) is mi = ni + 2n 2 + 1, m 2 = —n\ for the poles and 
mi = ni + 2n 2 , m 2 = —n\ for the zeros. 

Remark: let us consider the a\ — > oo limit of the one-parameter series R ai (see 
([35]). (1361) ') rewriting R ai (z) as R^iu) 

1 z 

Rb 1 (u) = R ai (z), with: 6 X = — , it = — . (83) 

Oi Oi 

In the ai — > oo limit, that is the bi — > limit, one easily verifies, order by order in 
u, that Rb 1 (u) becomes exactly the transcendental parametrization function (|48l) : 

Rbiiu) -> P(m) when &i 0. 

For ai = (—4)" (n — > oo), one finds that the radius of convergenc^ of the R ai (z) 
series becomes in the n — > oo limit R n ~ z s /4 n , in agreement with (|83|). 



K It is to the absolute value of the inverse of the image of by the n-th iterate of 3^, 4 of — 1. 
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2.5. P(z) and an infinite number of rational transformations: the sky is the limit. 

Note that some non- linear ODEs associated with P(z) and displayed in | Appendix B[ 
namely (|B.2[) and (|B.7[) . and the functional equation (1701 . are invariant by the 
change of variable (P(z), z) -> (-4P(z)/(l - P(z)) 2 , -Az). In fact pO)) . ([70]). 
and (|B.7|) are invariant by (P{z), z) -> (-4P(z)/(l - P(z)) 2 , -4z), but also 
(P(z), 0) -> (-(1 - P(z)) 2 /4/P(z), - z/4), and also by (P(z), z) -> (1/P(*), z). 
The function P(^) satisfies the functional equation: 

4P(z) , . 

pi-*-*) = -T-iw> ^ 

but also 

P((-7-24i)-z) = T(P(z)), (85) 
P((-7 + 24t).2) = T*(P(z)), 
and, more generally, as can be checked order by order on series expansions, 

P( ai -z) = P Q1 (p(z)). (86) 

For example, considering the "good" branch (|39p for the inverse of — 4z/(l — z) 2 , 
namely u(z), we can even check, order by order, on the series expansions of P(z) 

and S_^,^(z) the functional relation: 



S%(P(z)) = p(-|). (87) 



valid for \P(z)\ < 1 since the radius of convergence of S^^(z) is R = 1. 

Recalling the functional equations (|85p it is natural to say that if P(z) is singular 
at z = z 8 , then, for almost all the rational functions, in particular T(z) (resp. T*(z)) 
the T(P(z)) is also singular z — z s , and thus, from (J55J, P(z) is also singular at 
z = (— 7 ± 24 z) • z s . It is thus extremely natural to see the emergence of the infinite 
number of singularities in ([77jl of the form z — [N\ + i ■ AT 2 ) • z s , as a consequence 
of (|86|) together with a reduction of the one-parameter series R ai (z) to a rational 
function for an infinite number of selected values of a\ 1 namely the N\ + i ■ N2 in 
(|77[) . This is actually the case for all the values displayed in flTTJ) . For instance, for 
ai = 3 4 = 81 we get the following simple rational function: 

z 2 + 6z-3 \ 4 
v3 z 2 — 6 z — 1/ 

for which it is straightforward to verify that this rational transformation commutes 
with T(z), T*(z), —42/(1 — z) 2 , and is a solution of the Rota-Baxter-likc functional 
equation ([27)1 . The case ai = 5 4 = 625 in ([77)1 . is even simpler, since it just requires 
to compose T(z) and T*(z) 

R 6 2 5 (z) = T(T*(z)) = T*(T(z)) (89) 
/z 2 -2z + 5\ 4 /l - 12 z- 26 z 2 + 52 z 3 + z 4 \" 



^1 W = z -i o 2 V- 1 ) . ( 88 ) 



V5z 2 -2z+l/ Vl + 52z-26z 2 -12z 3 + z 4 

which, again verifies ()27[) and commutes with all the other rational functions, in 
particular (|88p. We also obtained the rational function corresponding to a\ = 7 4 = 
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2401, namely: 

iW*) = z-(^t\)\ with: (90) 

^2401 (A) y 

AW*) = z 12 -D 2 40i(~), and: (91) 

#2401(2) = 1 +196 z- 1302 z 2 + 14756 z 3 - 15673 z 4 - 42168 z 5 
+ 111916 z 6 - 82264 z 7 + 35231 z 8 - 19852 z 9 
+ 2954z 10 + 308z n -7z 12 . (92) 

The polynomial iV24oi(z) satisfies many functional equations, like, for instance 
(with R-i(z) = -4z/(l-z 2 )): 

4l2 '^ 2401 (iT^)) = d *m( z ) • D ™*{~) ( 93 ) 

and also: 

(l-z) 49 -^ 24 oi(i?-4W) 2 = D 240 i(z) 4 - z i9 ■ D 2i01 (^y . (94) 

We also obtained the rational function corresponding to a\ = ll 4 = 14641, 
namely: 

iWiW = ^( n 1464l( /1 ) 4 ' with; ( 95 ) 

,30 - 



N U6 4i(z) = * -#14641 - , and: (96) 



£»i464i (z) = 1 + 1210 z - 33033 z 2 + 2923492 z 3 + 5093605 z 4 



10 



385382514 z 5 + 3974726283 z 6 - 14323974808 z' 
57392757037 z 8 - 291359180310 z 9 + 948497199067 z 
1642552094436 z 11 + 1084042069649 z 12 + 1890240552750 z 13 
6610669151537 z 14 + 9712525647792 z 15 - 8608181312269 z ie 
5384207244702 z 17 - 3223489742187 z 18 + 2175830922716 z 19 
1197743580033 z 20 + 387221579866 z 21 - 50897017743 z 22 
7864445336 z 23 + 539 1 24 3 9 35 z 24 - 8 1 5 7 8 9 6 34 z 25 
28366041 z 26 - 50 9 2 9 56 z 27 + 207691 z 28 + 2 7 94 z 29 - 1 1 z 30 . 



and, of course, one can verify that i?i464i(z) actually commutes with i?_4, R&i, i?625, 
-R240i(z), and is a solution of the Rota-Baxter- like functional equation (|27[), Similarly 
to i?24oi(2) (see ([93]) . (|94p ). we also have the functional equations: 

4 30 - #14641 1 (z) ) = #14641 (*)• #14641 (~), (97) 



and also: 



(1 _ z) (4.30+i). ^ 14641 (i?_ 4 ( z )) 2 (98) 

L (z) 4 - z ( 4 - 39 +D.^ 14641 (I) 4 . 



= D 



146411 
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Next we obtained the rational function corresponding to a-i — 13 4 = 28561, 
which verifies (|2"T|) namely: 



R 



2856H 



{z) = z j*f™My } with: N2s56i(z) = z u. D!UB01 (l\ 



^28561 (2) = z i2 ■ ^28561 (J), and: (99) 

£>2856i(z) = (l-22z + 235z 2 - 228 2 3 + 39z 4 + 26z 5 + 13z 6 )-D^5 ) 61 (z) 

^28561 ( z ) = 1 + 2388 z- 61098 z 2 + 19225300 z 3 + 606593049 z 4 

- 1543922656 z 5 + 7856476560 z 6 - 221753896032 z 7 + 1621753072244 z s 

- 4542779886736 z 9 + 2731418674664 z 10 + 36717669656304 z 11 

- 200879613202428 z 12 + 547249607666784 z 13 - 934179604482832 z 14 

+ 1235038888776160 z 15 - 1788854212778642 z 16 + 3018407750933816 z 17 

- 4349780716415868 z 18 + 4419228090228152 z 19 - 2899766501472914 z 20 
+ 931940880451552 z 21 + 413258559018224 z 22 - 857795672629664 z 23 

+ 659989056851972 z 24 - 304241349909008 z 25 + 87636987790824 z 26 

- 14593362219920 z 27 + 1073204980340 z 28 + 45138167200 z 29 

- 23660433008 z 30 + 2028597792 z 31 - 29540327 z 32 + 32 3 8 4 20 z 33 

- 73386 z 34 - 492 z 35 + z 36 . (100) 

We get similar results, mutatis mutandis, than the ones previously obtained 
(commutation, functional equations like (|97l) . (|98)) ...). namely: 

4 42 - ^8561 (^T^y) - D 2S561 (z)-D 2S5ei {^j, ••• (101) 

The "palindromic" nature of (|88 ]l ([89 ]) . ([90 ]l . ([95 ]) and ([99 |) (see (|9T|) . (|96| ). (|9"9p 
corresponds to the fact that these rational transformations commute with J: 

1 



P , , ~oi(-), -^7^ - #625 (l) , ••• (102) 

In fact, more generally, we have Rfqi(l/z) = l/R N i(z) for N any odd integer 
(N = 9, 21, • • • , • • •) and Rn*{1/z) — Rj^i(z) for N any even integer. 

From (1771) one can reasonnably conjecture that the fourth power of any integer 
will provide a new example of R ai (z) being a rational function. The simple non-trivial 
example corresponds to the already found rational function: 

z- (1 - z) 2 

i? 16 (z) = 16- / J . (103) 
(z + 1) 

We already have explicit rational functions for all values of a\ of the form TV 4 for 
N = 2, 3, • • • , 16 and of course, we can in principle, build explicit rational functions 
for all the iV's product of the previous integers. Along this line it is worth noticing 
that the coefficients of the series R ai (z) are all integers when a\ is the fourth power 
of any integer. 

We are thus starting to build an infinite number of (elementary) commuting 
rational transformations, any composition of these (infinite number of) rational 
transformations giving rational transformations satisfying (|27l) and preserving the 
linear differential operator ft. This set of rational transformations is a pretty large 
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set! Actually this set of rational transformations corresponds to the isogenics of the 
underlying elliptic function. 

The proliferation of the singularities of P(z) corresponds to this (pretty large ...) 
set of rational transformations. Recalling (|87[) . the previous singularity argument is 
not vali^foi the (well-suited) inverse transformations ( S^, 4 (z), ...) of these rational 
transformations because (|57j) requires \P(z)\ < 1 (corresponding to the radius of 
convergence of S_y^(z)) and the singularity z = z s corresponds precisely to "hit" 
the value P(z) = 1. 



2.6. Other examples of selected Gauss hyper geometric ODE's 

For heuristic reasons we have focused on A(z) — (3 — 5 z)/z/(l — z)/4, but of course, 
one can find many other examples and try to generalize these examples. 
For instance, introducing 

1 rfln((l - z) 3 z 5 ) 1 5-8z 

(Z) ~ 6 Jz ~ 6 ' (l-z)z' 

the rational transformation 

R{z)= - 27 '(T^F' (104) 

verifies the "Rota-Baxter- like" functional relation (|2"T|) . This example corresponds to 
the following covariance [35] on a Gauss hypergeometric integral (of the c = 1 + b 
type, see below): 

2Fl ([il],[l}; z) = Z -^- f t-*"0-t)-V*.dt (105) 

-27- 



L 2'6 J L 6 J / 6 

' 1 -] [-] 

6 M 6 J ' "' (l-4z) 3 

which is associated with the elliptic curve: 

y 6 - {l-t) 3 -t 5 = 0. (106) 

Another example (of the c = 1 + a type, see below) is 

1 dln((l - zf z 2 ) 2 1- 2z , , 

3 «Z 3 (1 — Z) Z 

where the rational transformation 

R W = Z (\ ~. o ^3 = -8^- 36z 2 - 126z 3 - 387z 4 + ••• (108) 
(1 — 2 z) 

verifies the "Rota-Baxter- like" functional relation (|27|) . This example corresponds to 
the following covariance 25 on a Gauss hypergeometric integral: 

r , ,,1 2, r 4, \ z- 1 / 3 

2 -Pi 



1 2 - z /,1 2, r 4, z • (z - 2) 



2 1-22 V L 3' 3 J ' L 3 J ' (1 - 2z) 

| If this previous singularity argument were valid we would have had singularities as close as possible 
to 2 = (namely z s /(— 4) n ), yielding a zero radius of convergence. Similarly combining T*(z) and 
the inverse of T(z) we would have obtained an infinite number of singularities on the circle of radius 
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which is associated with the elliptic curve: 

y 3 - (1-i) 2 -t 2 = 0. (110) 

Note that, similarly to the main example of the paper, there exist many rational 
transformation:^ satisfying (f27|) that cannot be reduced to iterates of (|108[) . for 
instance: 

v3 



z ■ (1 - z) (z 2 - z + 1) 



T(z) = -27- ^ f- = -27 z -378z z (111) 

(z 3 + 3z 2 -6z + l) 3 

- 3888 z 3 - 34074 z 4 - 271620 z 5 - 2032209 z 6 + ■■■ 

One verifies immediately that (|111[) actually verifies (j2"7| with (|107|) . Not surprisingly, 
the two rational transformations (|108[) and (|lllj) commute. 

Another simple example with rational symmetries corresponds to — (D z + 
A(z)) ■ D z with 

w 2 z(l-z) 2 dz v 7 

It has the simple (genus zero) hypergeometric solutior jftl : 

F{z) = z 1 / 2 ■ a Fi([l, d'tf];*) = arcian/i(z 1 / 2 ). (113) 

The linear differential operator f2 is covariant under the change of variable z — > l/z 
and z — > i?(z) wher^J: 

i?(z) = --4^- (H4) 
(1 + z)^ 

One can easily check that (|1 12[) and (| 1 14[) satisfy the functional equation (|27|) . One 
also verifies that (|112l) and z — >• l/z or the iterates of (II 14[) satisfy the functional 
equation ([27]) . The solution of the adjoint operator are (1 — z) • z 1 / 2 and 

F(z) = z-(l-z). 2 Jl([l, (H5) 

2 2 

= z 1/2 • (1 - z) • arctanh(z 1/2 ) = z z 2 z 3 

3 15 

-—z 4 -—z 5 -—z 6 -—z 7 + ■■■ 
35 63 99 143 + 

One verifies, again, that (|115[) and (| 1 14[) commute, (|1 15|) corresponding to the 

"infinitesimal composition" of (|1 14|) (see (|5T|)). 

A first natural generalization amounts to keeping the remarkable factorization 
(|16[) which will, in fact, reduce the covariance of a second order operator to the 
covariance of a first order operatoiSJ Such a situation occurs for Gauss hypergeometric 
functions 2Fi([a, 6], [1 + a]; z^ solution of the (a, &)-symmetric linear differential 
operator 

z-(l-z)-!? 2 + (c- (a + b+1) ■ z) ■ D z - a- b, (116) 

| Note a (small) misprint in formula (64) page 174 of Vidunas 25 . 
ff Of the c = 1 + b type (see below). 

f The change of variable II114H can be parametrized with hyperbolic tangents: z — > z' with 
z = tanh(u) 2 , z' = tanh(2w) 2 . Note that z — > 4-z/(l — z) 2 is parametrized by z = tan(u) 2 , z' = 
tan(2u) 2 but z — ► —4 ■ z/(l — z) 2 is not parametrized by trigonometric functions. 
% Thus avoiding the full complexity (and subtleties) of the covariance of ODE's by algebraic 
transformations like modular transformations J51 . 
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as soon aqj c = 1 + a. For instance 

F(z) = z a ■ 2 F 1 ([a, b], [1 +o]; zj, (117) 

is an integral of a simple algebraic function and is solution with the constant function 
of the second order operator: 

/ (a — b — 1) z + 1 — a\ , 
°=(^ + z.(l-z) H < 118 > 

= (d, + d ^-?-*-Kp» 

V dz J 

yielding a new A{z): 

= (l- a )+( a -b-l)z = l_a 

(1 — z) ■ Z Z 1—2 

The adjoint of (111 8[) has the simple solution z 1 ~ a ■ (1 — z) h : 

F{z)= z-[l-z) b - 2 F 1 ([a, b], [1 +o]; zj. (120) 

Due to the (a, &)-symmetry of ()116|) we have a similar result for c = 1+6. The 
function J-{z) = z b ■ 2-F1 ([a, b], [1 +b]; z^j is solution of (|118j) where a and b have 
been permuted: 

(b- a- 1) z + 1 - b\ , 
z • (1 — z) / 
yielding another 

(1 - z) • z 

The adjoint of (|121j) has the solution (1 — z) a ■ z 1 ^ together with the hypergeometric 
function: 

F(z) = z-(l-z) a - 2 F 1 ([a,b],[l + b]; *). (123) 

where one recovers the previous result (| 1 15|) . 

We are seeking for (Gauss hypergeometric) second order differential equations^) 
with an infinite number of (hopefully rational, if not algebraic) symmetries: this 
is another way to say that we are not looking for generic Gauss hypergeometric 
differential equations, but Gauss hypergeometric differential equations related to 
elliptic curves, and thus having an infinite set of such isogenics. We are necessarily 
in the framework where the two parameters a and b of the Gauss hypergeometric 
are rational numbers in order to have integral of algebraic functions (yielding globally 
nilpotent [12] second order differential operators). Let us denote by D the common 
denominator of the two rational numbers a — N a /D and b — N^/D, the function 
(|117[) is associated to a period of the algebraic curve: 

y D = (l-t) N *-t D - N °. (124) 

X See for instance iC.2i in |Appcndix C| 

§ More generally in our models of lattice statistical mechanics (or enumerative combinatorics etc.) 
we are seeking for (high order) globally nilpotent 1221 operators that, in fact, factor into globally 
nilpotent operators of smaller order, which, for Yang-Baxter integrable models with a canonical 
elliptic parametrization, must necessarily "be associated with elliptic curves". | Appendix D| provides 
some calculations showing that the integral for x' 2 ^ the two-particle contribution of the susceptibility 
of the Ising model 31 . 32 , 33] is clearly, and straightforwardly, associated with an elliptic curve. 



Renormalization, Hauptmoduls 



26 



-1/6 



t~ 5 / 6 (1 - t)~V3 . dt, 



We just need to restrict to triplets of integers (N a , Nf,, D) such that the previous 
curve is an elliptic curve. 

Let us give an example (of the c = 1 + b type) that does not correspond to a 
genus one curve, with 

[I h rli ^ I 

[ 3' 6 M 6 J ' J ~ 6 
which corresponds to the genus two curve: 

y 6 - (l-t) 2 -t 5 = 0. (125) 
Again one introduces A(z): 

1 dln((l - z) 2 z 5 ) 1 5-7 z 
(Z) ~ 6 Tz ~ 6 ' z-(l-z)' 

and seek for R(z) as series expansions analytical at z = 0. One gets actually, order 
by order, a one-parameter family: 

2 

R ai (z) — a\ ■ z — — oi ■ (ai — 1) • z (126) 
+ g^ai .(oi-l).(17ai -87)-z 3 

' 2 or. , ,<oo\ „4 



■ ai • (ai - 1) • (113 a{ - 856 a x + 3438) • z 



84721 

■ oi • foi - 1) • (3674 a? + 121194 ai 
38548055 v ' v 1 1 

- 552261 ai + 2095059) ■ z 5 + ■■■ 

+ l^- a i-( a i -l)--Pn(ai)-^ + •■■ 

where e„ = for n odd and e„ = 1 for n even, and N(n) is a (large) integer 
depending on n, and P„(ai) is a polynomial with integer coefficients of degree 
n — 2. One easily verifies, order by order, that one gets a one-parameter family 
of transformations commuting for different values of the parameter: 



R, 



ai (R bl (zfj = R bl (^R ai (z)j = R aibl (z). 



As far as the "algorithmic complexity" of this series (|126|) is concerned it is worth 
noticing that the degree growth |M] of the series coefficients, is actually linear and not 
exponential as we could expect [35 at first sight. Even if this series is transcendental, 
it is not a "wild" series. 

Seeking for selected values of ai such that the previous series (|126l) reduces 
to a rational function one can try to reproduce the simple calculations (14"4"|) . (05]), 
but unfortunately "shooting in the dark" because we have no hint of a well-suited 
denominator (if any!) like the polynomials in the lhs of (|44|) . (j45|) . 

It is also worth noticing that if we slightly change A(z) into: 

the algebraic curve (|125[) becomes y N — (1 — t) 2 ■ t 5 — which has, for instance 
genus five for N = 11, but genus zero for TV = 7. For any of these cases of (|127l) 
one can easily get, order by order, a one-parameter series R ai totally similar to (|126p 
with, again, polynomials P n (a\) of degree n — 2. 
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The first coefficient a 2 is, in general: 

«2 = -^TT5" -!)■ (128) 

For the genus zero case, N = 7: 

a 2 = ~ ■ a x ■ (a x - 1), a 3 = - jjjjg • «i • («i - 1) ' (127 - ax), 

a+ = -T^f^' a i-( a i -l)'(254a? + 185 oi +7499), ••• 
ax ■ (a x - 1) ■ P n (ax). 



N(n 

which corresponds to the solution: 



y J\-^.{l-z)-^.dt = z 2 ' 7 - ^([2,2], [2] 



(129) 



Using the parametrization of the genus zero curve 

(u + l) 2 -u 5 _ u 7 

y ~ ~{u + iy - u 7 ' 1 ~ ~ (u + iy - u 7 ' 

one can actually perform the integration (II 29[) of dt/y and get an alternative form of 
the hypergeometric function (|129[) : 

z~ 5/7 ■ (1 - z)- 2/7 ■ dt = / p(u)-du = / =-dv, 

Jo Jo 1 - v 7 

. u 7 (u + l) 4 • u 

where: z = —-, -r= ~, p(u) — 



(u + 1) 7 - u 7 ' ry ' (u+ 1) 7 

and: v = 



7 

it IT 



1 + u w 7 -l 
Except transformations like v — > u • v (with oj 7 = 1) which have no impact on z, it 
seems difficult to find rational symmetries in this genus zero case. 
For N = 11 (genus five) the first successive coefficients read: 
2 

a 2 = -— ■ ax ■ (ax - 1), 

° 3 = ~80^2 ' Ql ' ~ ^ ' ( 143ai + 367 )' 

fl 4 = -7^7^- a i-( a i -!)• (H86a? + 2473a! +5011), ••• 
206346 

■ax - (ax - 1) ■ P n (ax). 



N(n) 

The "infinitesimal composition" function F(z) (see (|5Tj) . (|52|) and (f53| ) reads: 

, . <9i? ai , 2 2 15 3 5 4 37 5 

F(z) = — - — -ai-l = z Z z z z 

dax 11 17 238 119 1190 

6 2183 7 4366 8 58941 9 



(130) 



36295 108885 258213 4045337 
1807524 ^ 46543743 ^ n 5305986702 



141586795 4106017055 521464165985 

and, again we can actually check that this is actually the series expansion of the 
hypergeometric function 

z-(l-z)- 2 Fx([l, ^j], (131) 
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solution of 57* the adjoint of the f2 linear differential operator corresponding to this 
(genus 5) N — 11 case: 



We have similar results for (|lf 7[) . (|118p . (|119l) . As far as these one-parameter 
family of transformations R ai , are concerned, the only difference between the generic 
cases corresponding to arbitrary genus and genus one cases like (|107p . is that, in the 
generic higher genus case, only a finite number of values of the parameter a\ can 
correspond to rational functions. Note that this higher genus result generalizes to the 
arbitrary genus Gauss hypergeometric functions (| 1 1 7[) and associated operators (|118|) 
and function (|119p . In this general case one can also get order by order a one-parameter 
family of transformations R ai satisfying a commutation relation (11271) . 

Note that R(z) = l/z is actually a solution of (|2T[) for this genus-two example 
(I126p . Along this line of selected R(z) solutions of (|27|) many interesting subcases of 
this general case (|117p . (|118p . (|119p are given in | Appendix C.l| 

In our previous genus-one examples, with this close identification between the 
renormalization group and the isogenies of elliptic curves, we saw that, in order to 
obtain linear differential operators covariant by an infinite number of transformations 
(rational or algebraic), we must restrict our second order Gauss hypergeometric 
differential operator to Gauss hypergeometric associated to elliptic curves (see 
|Appcndix C| and |Appendix Dp . Beyond this framework we still have one-parameter 
families (see (I127p ) but we cannot expect an infinite number of rational (and probably 
algebraic) transformations to be particular cases of such families of transcendental 
transformations. 

3. Conclusion 

We have shown that several selected Gauss hypergeometric linear differential operators 
associated to elliptic curves and factorised into order one differential operators 
actually present an infinite number of rational symmetries that actually identify 
with the isogenies of the associated elliptic curves that are perfect illustrations 
of exact representations of the renormalization group. We actually displayed all 
these calculations, results and structures because they are perfect examples of exact 
renormalization transformations. For more realistic models (corresponding to Yang- 
Baxter models with elliptic parametrizations) , the previous calculations and structures 
become more involved and subtle, the previous rational transformations being replaced 
by algebraic transformations corresponding to modular curves. For instance, in 
our models of lattice statistical mechanics (or enumerative combinatorics, etc.), we 
are often getting globally nilpotent linear differential operators [22] of quite high 
orders [351 l3"7]_[3"8l I3U1 HP] that, in fact, factor into globally nilpotent operators of 
smaller ordersjfl, which, for Yang-Baxter integrable models with a canonical elliptic 
parametrization, must necessarily "be associated with elliptic curves". |Appendix D| 
provides some calculations showing that the integral for X ■> the two-particle 
contribution of the susceptibility of the Ising model [311 [32l [33] is clearly and 
straightforwardly associated with an elliptic curve. 

K Experimentally 1371 and as could be expected from Dwork's conjecture 1221 . one often finds for these 
small order factors hypergeometric second order operators and sometimes selected Heun functions 141 1 
(or their symmetric products). 




(132) 



Renormalization, Hauptmoduls 



29 



We wanted to highlight the importance of explicit constructions in answering 
difficult or subtle questions. 

All the calculations displayed in this paper are elementary calculations given 
explicitly for heuristic reasons. The simple calculations (in particular with the 
introduction of a simple Rota-Baxter like functional equation) should be seen as some 
undergraduate training to more realistic renormalization calculations that will require 
a serious knowledge of fundamental modular curves, modular forms, Hauptmoduls, 
Gauss-Manin or Picard-Fuchs structures [42l [43] and, beyond, some knowledge 
of mirror symmetries [T51 [T7J [TSJ [T|3 [2D] of Calabi-Yau manifolds, these mirror 
symmetries generalizing^ the Hauptmodul structure for elliptic curves. 
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Appendix A. Comment on the Rota-Baxter-like functional equation (|27|) 



We saw, several times, that the Rota-Baxter-like functional equation ([27]) is such that 
for a given A(z) one gets a one-parameter family of analytical functions R(z) obtained 
order by order by series expansion (see (|36p . (11261) ). Conversely for a given R(z), for 
instance R(z) = — 4z/(l — z) 2 , let us see if R(z) can come from a unique A(z). 
Assume that there are two A(z) satisfying (|30|) with the same R(z) = — 4 z/(l — z) 2 . 
We will denote S(z) the difference of these two A(z), and we will also introduce 
A(z) — z ■ 8(z). It is a straightforward calculation to see that A{z) verifies 

. , ^ 1 + Z . / —4 Z \ f . . 

= i-r A ((i^F)< ( ) 

which has, beyond A(z) = 0, at least one solution analytical at z — that we can 
get order by order: 

. , , 2 22 2 394 , 262634 . 

A(z) = 1 + - z H z 2 H z 3 H z 4 + ■■■ 

w 5 75 1625 1243125 

It is straightforward to show from (|A.1[) . from similar arguments we introduced for 

(|52p on the functional equations (|53"|) that A(z) is a transcendental function. 

Appendix B. Miscellaneous non-linear ODE's on P(z) 
From (p)2"j) one can get 

P(2) 

F'(P(z)) = 1 + z-p^, 

§ See for instance equation (1.9) of |17 | . Do note that the periods of certain K3 families (and 
hence the original Calabi-Yau family) can be described by the squares of the periods of the elliptic 
curves |17| . The mirror maps of some K3 surface families are always reciprocals of some McKay- 
Thompson series associated to the Monstruous Moonshine list of Conway and Norton, the mirror 
maps of these examples being always automorphic functions for genus zero |18l I19 |. 
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p(2) p(3) (P(2h2 

F" (PCs 11 = — 4- z ■ — - z ■ — — 

V\ Z >> (p(l))2 +Z (P(l))2 Z (p(D)3' 

and from ([57)) . the linear second order ODE on F(z), one deduces the third order 
non-linear ODe(3 on the (at first sight non-holonomic) function P(z): 

z • (5 P 2 - 6 P + 3) • (P (1) ) 4 - P ■ (5 P - 3) • (P - 1) ■ (P (1) ) 3 
- z ■ (P - 1) • P • (5P - 3) ■ P (2) ■ (P (1) ) 2 
+ 4P 2 • (P - l) 2 ■ (P< 2 ) + 2 ■ P( 3 )) • P« (B.l) 
-4z- (P< 2 )) 2 • P 2 • (P- l) 2 = 0, 

where the P' n )'s denote the n-th derivative of P(z). This third order non-linear ODE 
has a rcscaling symmetry z — > p ■ z, for any p, and, also, an interesting symmetry, 
namely an invariance by z — > z a , for anjjj] value of a. 

In a second step, using differential algebra tools, and, more specifically, the fact 
that P(Q(z)) — Q(P(z)) — z together with the linear ODE for Q(z), one finds the 
simpler second order non-linear ODE for P(z): 

^ (2) -r^iyp-( p(1) ) 2 + rr p(1) = °« (R2) 

or 

P (2) _ ( 3 . 1 + 1 . _J_^ . ( P (i))2 + 3.1. p(i) = o. 
U P 2 P - 1/ v j 4 z 

Note that, more generally, the second order non-linear ODE 

P< 2 > -fi.l + i. -L-) . (p(D)2 + 1 . p(D = o, (B.3) 

V4 P 2 P - 1/ V ' z ' v; 

yields (|B.1[) for any value of the constant 77. The change of variable z ^ z™, changes 
the parameter 77 into 1 + a • (77 — 1). In particular the involution z <-> I/2 changes 
7/ = 3/4 into 77 = 5/4. 

This non-linear ODE, looking like Painleve V, is actually invariant by the change 
of variable P — > — 4P/(1 — P) 2 . It is, also, invariant by any rescaling z -> Az, like 
the particular degenerated subcase of Painleve V: 

y» - (±- +^—).y' 2 +1V = 0. (B.4) 

\2y y — 1/ z 

With (|f)2")) we recover the "Gauss-Manin" idea of Painleve functions being seen 
as deformations of elliptic functions: 

Z -^T = Pi*) ■ 0- Pi*)) 1 * • lu^;P(z))- (B.5) 

or: 

d arctanh((l - P(z)) 1 / 2 ) -, 1 r 

hu "= ^([ 3 . 5 Mi]^W 

In fact, recalling Q(P(z)) = z, one also has the relation 

P(z). 2 P 1 ([i i],[^] ; p(z)) 4 = z, (B, 

f Using differential algebra tools one can verify that l|75j l implies I)B,1J| . 
f Beyond diffeomorphisms of the circle: the parameter a can be a complex number. 
§ Having the movable-poles solutions: (a@ +z^) 2 /(a^ —z@) 2 . 
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yielding with (|B.5|) the simple non-linear order-one differential equation 



~ 3 (P') 4 -(l-P) 2 -P 3 = 0, (B.7) 



z 



already seen with (1741) . and that we can write in a separate way: 
dP dz 



(B.8) 



(1 - py/2 . p3/4 z 3/4- 

Note that P{z i - { - 1 -^) is actually solution of (|B~3|) . 

Equation (|B.7[) has (|B.6|) as a solution wbut in general the Puiseux series solutions 
Pa{z) of the functional equation (/i is a constant): 

P,{z) 1/A - 2F 1 {\±,\u\]\ Pa{z)) = M +^ 1/4 , 
P,(z) = Pfa + z 1 '*) 4 ), 

P^z) = /V)+ 4 • ^ 3 ■ PV) ' ^ 1/4 + ••• (B.9) 

It is a straightforward exercise of differential algebra to see that the order-one non- 
linear differential equation (|B.7[) implies (|B.2[) . In particular not only (|B.6|) is solution 
of (|B.2[) but also all the Puiseux series solutions (|B.9[) of (|B.7[) . More generally the 
solutions of the functional equation: 

P*\(z) 1/A - 2F 1 {[\, PkxW) = A^+A-z 1 ^. (B . 10) 

verify (IB.2|) . This corresponds to the fact that 

z 3 -(P') 4 - A 4 -(1- Pf ■ P 3 = 0, (B.ll) 
yields (|B.1|) which is scaling symmetric (z —> p- z) when (|B.7[) is not. More generally 

z iv.r P >y -A 4 -(l-P) 2 -P 3 = 0, (B.12) 

yields (|B.1[) for any value of the parameters n and A. Finally, one also has that the 
solution of the functional equation 

P^z) 1 ^ ■ 2 F 1 ([\,\],[^ P v {z)) = a^A-z 1 -", (B.13) 

are solutions of (jB.ll) . but also of (|B.3I) and even of ()B.12|) . 

Equation (|B.3[) with rj = 1/2 (instead of 77 = 3/2 in (|B.2[) ). has a solution, 
analytical at z = : 

1 2 7 3 1 4 121 5 7 6 
1 + x -\ — x H a; H x -\ x ti -\ x b 

2 40 20 9600 2400 

211 7 41 o , . 

+ M S+ 3i20i I + (R14) 
This series has a singularity at — 1/4 • z"^ where z s is given by (|76|) . The radius of 
convergence of (|B. 14|) corresponds to this singularity, namely R — 1/4 • z 2 . This 
singularity result can be understood from the fact , at r\ = 1/2, P(z 2 ) is actually 
solution of (|B.3|) . 

In fact, we have the following solutions of (|B.3|) for various selected values of rj. 
For r\ = 0, P(z 4 ) is solution of (lEO]) . For 77 = 2/3, P(z 4 / 3 ) is solution of (|B~3|) . 
and, more generally, P(z 4 '( 1_,? ^) is solution of (|B.3j) . 



Renormalization, Hauptmoduls 32 
Appendix C. Gauss hypergeometric ODE's related to elliptic curves 



It is not necessary to recall the close connection between Gauss hypergeometric 
functions and elliptic curves, or even modular curves [HI |3S] and Hauptmoduls. This 
is very clear on the Goursat-type relation 

(, 2a+ 1, r 4a + 2, \ , s 

2 F 1 ([2a,— ] ;a: ) = (c.l) 

2 / a a + 1 r 4a + 5 1 27 (g - l) 2 ■ x 2 x 

which generalizes the simpler quadratic Gauss relation: 

_ /. ,a + b+l \ f a b a + b+l \ 
2 F l ([a,b],[ };xj = 2 F 1 {[-,-],[ ];4a; (1 - x)j. 

On (jC.ip one recognizes (the inverse of) the Klein modular invariant for the pull-back 
of the hypergeometric function on the rhs. 

Many values of [[a, b], [c]] are known to correspond to elliptic curves like 
[[1/2, 1/2], [1]] (complete elliptic integrals of the first and second kind) or modular 
forms: [[1/12, 5/12], [1]], [[2/3, 2/3], [1]], [[2/3, 2/3], [3/2]] and they can even be simply 
related: 

'7 4- 27\ V 3 

" 1 • 2 *i([2/3,2/3],[l]; -1/27*) = 



M (z) • 2 *i([l/12, 5/12], [1]; 1728 * 

v (z + 27)(z + 3) 

(z + 27) (z + 3) 3 \-Vi2 



27 



where: u(z) 

PV ' V 729 

Once we have a hypergeometric function corresponding to an elliptic curve for 

some values of (a, 6, c), one can find other values of (a, b, c) also corresponding to 

elliptic curves: 

a fi([o,6],[c];x) — > .x 1 ^ • 2 F X ([1 + a - c, 1 + 6 - c], [2 - c]; 

In order to provide simple examples of linear differential ODE's we will restrict 
ourselves (just for heuristic reasons) to Gauss hypergeometric second order differential 
equations. 

Let us recall the Euler integral representation of the Gauss hypergeometric 
functions: 

2 F 1 ([a, b], [c]; z) (C.2) 



r(c) 




Z" 1 du> 




(1- 


-w) c -^ b • (1- 


r(6)r( c - 


6) 


Vo w 


r(c) 






w a 


•(1 


-wf- 1 - 1 • (1 


T(a)T(c- 


a) 


io w 


r(c) 




Jo 


du 


u a • 


(1 _«)c-l-a 


i»r(c- 


a) 


u 



v-6 



(1-u)- 



On the last line of (|C.2[) . the selected role of c = 1 + a is quite clear. 

f Taking for x the elliptic lambda function. 
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Recall that the corresponding second order differential operator is invariant 
under the permutation of a and b which is not obvious^ on the Euler integral 
representations of the hypergeometric functions (this amounts to permuting and 
oo). The permutation of a and 6 is always floating around in this paper. 

When the three parameters a, b and c of the Gauss hypergeometric functions 
are rational numbers we have integrals of algebraic functions and, therefore, we 
know [22l SSI |47l 0H1 [49] that the corresponding second order differential operator 
is necessarily globally nilpotent [22j [46] [4JJ |48l [49]. Let us restrict to a, b and c being 
rational numbers a — N a /D, b = N^/D and c = N c /D where D is the common 
denominator of these three rational numbers. The Gauss hypergeometric functions 
are naturally associated to the pencil of algebraic curves: 

y°= (l-u) N >.u°-"».(l-^y Nc+D+N \ (C.3) 



Recalling the main example of the paper, one associates with 2-F1 ( [1/4, 1/2], [5/4]; z 

2^1 ([1/2, 1/4], [5/4]; z) = 
r(5/4) f 1 dw 



r(i/2)r(3/4) 
r(5/4) 



w 1 ' 2 ■ (1 - ™r 1/4 • (1 - zw)- 1 ^ 
- 1 ' 2 ■ f u- 1 ' 2 ■ (1 - -r 1 / 4 • (1 - u)- 1 / 4 ■ du 

Jo z 



r(i/2)r(3/4) 

the z-pencil of elliptic curvetQ 

y 4 - u 2 ■ (1 - u) ■ (1 - -) - 0, (C.4) 

z 

where we associated (see (fT2T ) to 2 -Pi([l/2, 1/4], [5/4]; zj the elliptic curve 

y 4 - u 3 ■ (1 - u) 2 = 0. (C.5) 



Appendix C.l. Miscellaneous examples 

In the more general (fTTTj) . (|TT5|) . pig|) . (resp. ([T2T]) , {T25])) framework, one can find 
many interesting subcases. 

• The previous R(z) = 1/z involution is solution of the functional relation (|27|) 
when a = 2 6 if c = 1 + b, or b = 2a if c = 1 + a. 

• The involution i?(z) = 1 — 2 is solution of the functional relation ([27]) when 
a + &= lif c = 1 + 6, or c = 1+a. 

• The infinite order transformation: 

R(z) = t — , R {n \z) = t n r , 

w 1 + (t- l)-z' w 1 + (t n - l)-z 

is solution of the functional relation (|27j) when a = 1 + 6 if c = 1+6, or b = 1 + a 
c = 1+a. 

• The scaling transformation R(z) = t ■ z is solution of the functional relation 
(1271) when a = and c = 1 + 6 (resp. 6=0 and c = 1 + a). 

I For instance for 2Fi([l/4, 1/2], [5/4]; 2) it changes an Euler integral with r(i/4(r^(i) = 3 ^ n *° 
an Euler integral with ^^^L^ = J . . 

f The algebraic curves IIC.4I I are genus one curves for any value of 2, except 2=1 where the curve 
becomes the union of two rational curves (u 2 — u + y 2 ) (u 2 — u — y 2 ) = 0. 
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• We also have a quite degenerate situation for b = 1 or a = 1 when c = 2 
with the infinite order transformation 

R(z) = 1 -t-(l-z), R {n) (z) = 1 -t n -(l-z), 

solution of ([2"7]) . 

• The two order-three transformations 

R{z) = ^— !-, R{R{z)) = ' 



z 1 — z 

are solutions of the functional relation (f2"T| for a = 2/3, = 1/3, c = 4/3, or 
a = 1/3, b = 1/3, c = 4/3. 

Appendix D. Ising model susceptibility : x^ and elliptic curves 

The two-particle contribution of the susceptibility of the Ising model [3TJ [32[ [33J is 
given by a double integral. This double integral on two angles x*- 2 ) reduces to a simple 
integral (because the two angles are opposite): 

_ (2) = r i+* 2 f*-M*)\> 

Jo 

where 



1 — x 2 V 1 — x 



x = A — B, A = J- - cos(f), B 2 = A 2 - 1. 

2 w 

2 1 

2/ = 



A 2 - I' 

Denoting C = cos(O) we can rewrite the integral x^ as : 

~(2) _ f 1 dC 2 2 1 + x 2 (r> 11 

x J (i-c 2 y/ 2 ' x ' v ' (i -x 2 ) 3 ' 1 j 

that we want to see as: 

" ^ ^ * (D.2) 



z Jo Z 

The variable z reads: 

1 1 2 2 1 + 2:2 n , , 

- ~ (1 _ C 2)i/2 - x ^ ■ (1 _ X 2 )3 - ( D - 3 ) 

which after simplifications gives 

A 2 (C* 2 - 1) ■ z 2 + (A 2 - l) 5 = 0, (D.4) 

that is 

(± -Q 2 -(C 2 -l)-z 2 + ((-L -C) 2 -If = 0. 
2w 2w 

In terms of the variable q — w ■ C one can rewrite (see (|D.2|) ^ the integral (|D. 1|) as 

an incomplete integral: 

256 -(1- 2 q) 2 (q 2 -w 2 )- Z 2 w A 

+ (2q-l + 2w) 5 (2q-l-2wf = 0. 

This w-pencil of algebraic curves is actually a w-pencil of genus one curves, seen as 
algebraic curves in Z and q. 

| The prefactors in front of the integrals are not relevant for our discussion here. 
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